RENORMALIZATION GROUP IN THE STATISTICAL THEORY OF 
TURBULENCE: TWO-LOOP APPROXIMATION 
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The field theoretic renormalization group is applied to the stochastic Navier-Stokes equation that 
describes fully developed fluid turbulence. The complete two-loop calculation of the renormalization 
constant, the beta function and the fixed point is performed. The ultraviolet correction exponent, 
the Kolmogorov constant and the inertial-range skewness factor are derived to second order of the 
e expansion. Submitted to Acta Physica Slovaca. 
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1 I. INTRODUCTION 

One of the oldest open problems in theoretical physics is that of describing fully developed turbulence on the basis 
of a microscopic model. The latter is usually taken to be the stochastic Navier-Stokes (NS) equation subject to an 
external random force which mimics the energy input by the large-scale modes; see, e.g., jl],|j. The aim of the theory 
is to verify the basic principles of the celebrated Kolmogorov-Obukhov phenomenological theory, study deviations 
^""^ ' from this theory, determine the dependence of various correlation functions on the times, distances, external (integral) 
and internal (viscous) turbulence scales, and derive the corresponding scaling dimensions. Most results of this kind 
were obtained within the framework of numerous semiphenomenological models which cannot be considered to be the 
basis for construction of a regular expansion in certain small (at least formal) parameter ^,|^| . 

An important exception is provided by the renormalization group (RG) method that was earlier successfully applied 
in the theory of critical behaviour to explain the origin of critical scaling and to calculate universal quantities (critical 
dimensions and scaling functions) in the form of the e expansions || . 

The RG approach to the stochastic NS equation, pioneered in [4-7] , allows one to prove the existence of the infrared 
(IR) scale invariance with exactly known "Kolmogorov" dimensions and the independence of the correlation functions 
of the viscous scale (the second Kolmogorov hypothesis), and calculate a number of representative constants in a 
reasonable agreement with experiment. Detailed review of the RG theory of turbulence and more references can be 
found in gg]. 

In contrast to the standard 4 model of critical behaviour , where the critical exponents are known up to the order 
{SJ . e: 5 (five-loop approximation), all the calculations in the RG approach to the stochastic NS equation have been confined 
with the simplest one-loop approximation. The reason for this distinction is twofold. First, the multiloop calculations 
for this dynamical model are rather involved: one can say that the two-loop calculation for the stochastic NS equation 
is as cumbersome as the four-loop calculation for the conventional </> 4 model. Second, the critical dimensions for the 
most important physical quantities (velocity and its powers, frequency, energy dissipation rate and so on) are given by 
the one-loop approximation exactly (the corresponding e series terminate at first-order terms) and the higher-order 
calculations for them are not needed. 

However, the e series for other important quantities do not terminate and the calculation of the higher-order terms 
for them is of great interest. In this paper, we present the results of the two-loop calculation for a number of such 
quantities: the ultraviolet (UV) correction exponent ui, the Kolmogorov constant Ck and the inertial-range skewness 
factor S. 



II. THE MODEL, FIELD THEORETIC FORMULATION AND RENORMALIZATION 

Detailed exposition of the RG theory of turbulence and the bibliography can be found in below we restrict 

ourselves to only the necessary information. 

As the microscopic model of the fully developed, homogeneous, isotropic turbulence of an incompressible viscous 
fluid one usually takes the stochastic NS equation with a random driving force 

Vm =v d 2 <pi -diP + Fi, V t =d t + { V d). (2.1) 

Here <pi is the transverse (due to the incompressibility) three-dimensional vector velocity field, P and Fi are the 
pressure and the transverse random force per unit mass (all these quantities depend oniE {t, x}), uq is the kinematical 
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viscosity coefficient, d 2 is the Lapface operator and Vt is the Lagrangian derivative. The problem (2.1) is studied on 
the entire t axis and is augmented by the retardation condition and the condition that ipi vanishes for t — > — oo. We 
assume for F a Gaussian distribution with zero mean and correlator 

(Fi^Fjix')) =5(t-t')(2Tr)- 3 J dkP y (k)d F (fc) exp [ik(x - x')] , (2.2) 

where -Pij(k) = Oy — kjkj/ k 2 is t he t ransverse projector and dp(k) is some function of k = |k| and model parameters. 

The stochastic problem (2.1), (2.2) is equivalent to the field theoretic model of the doubled set of fields $ = {ip, ip'} 
with action functional 

S($) = p'D F p>'/2 + <f/[-dt<p + v Q d 2 V - (ipd)<p], (2.3) 

where Dp is the random force correlator ( [2.2j ) and the required integrations over x = {£, x} and summations over the 
vector indices are understood. 

The standard RG technique can be applied to the model (^J^) if the function dp(k) is chosen in the form 

d F (k) = Dok 1 - 2 " h(m/k), h(0) = l. (2.4) 

Here Do is the positive amplitude factor, m — l/L is the reciprocal of the integral turbulence scale L, the function 
h(m/k) provides the IR regularization and the exponent e > plays the part of the RG expansion parameter, similar 
to that played by e = 4 — d in Wilson's theory of critical phenomena . The real (physical) value of this parameter 



is e = 2: idealized energy injection by infinitely large eddies corresponds to dp{k) oc <5(k), and the function (2.4) for 
e — > 2 and the appropriate choice of the amplitude can be considered as a power-law model of the three-dimensional 
6 function. 

The model (2.3) is logarithmic (the coupling constant go = Dq/v^ is dimensionless) at e — 0, and the UV divergences 



have the form of the poles in e in the correlation functions of the fields ip and tp' . Superficial UV divergences, 
whose removal requires counterterms, are present only in the 1-irrcducible function {p'p), and the corresponding 
counterterm has the form ip'd 2 p. Thus for the complete elimination of the UV divergences it is sufficient to perform 
the multiplicative renormalization of the parameters Uq and go — Dq/vo with the only independent renormalization 
constant Z v : 

u = uZ v , g = g f i 2£ Z g , Z g = Z~ 3 {D = g u 3 = gp? e v 3 ). (2.5) 

Here \± is the reference mass in the minimal subtraction (MS) scheme, which we always use in what follows, g and v 
are renormalized analogues of the bare parameters go and i>o, and Z = Z(g,£,d) are the renormalization constants. 
In the MS scheme they have the form "1 + only poles in e," in particular, 

oo oo n 

Z v = 1 + a k (g)£- k = 1 + 9 n a ^ £ " fe ' ( 2 - 6 ) 

fc=l n=l k=l 

with the one- loop coefficient an = —l/40n 2 

III. TWO-LOOP APPROXIMATION FOR THE RG FUNCTIONS, FIXED POINT AND THE UV 

CORRECTION EXPONENT 

We have performed calculation of the constant Z v with the accuracy of 0(g 2 ) (two-loop approximation). The 
calculation is rather involved and will be presented elsewhere (some details can be found in ju|), and below we giv e 
only the results for the residues a 2 2 and a 2 i at the second-order and first-order poles in e in the representation (p^q): 

a22/a\i = 1, a 2 i/al 1 ~ -1.65. (3.1) 

The knowledge of the renormalization constant Z v to order 0{g 2 ) allows for the calculation of the RG functions, 
the anomalous dimension ^ v and the beta function f3 g (£,g), with the following accuracy: 

0(g,e)=V li g = g(-2E + &y v (g)), 
7„($) = T>» \nZ v = -2gd g ai(g) = -2 (a n g + 2a 21 g 2 ) + 0(g 3 ), (3.2) 



2 



where is the operation /id/d/i at fixed bare parameters. In the MS scheme only the residues at the first-order 
poles in e, that is, only the coefficients aki, contribute to the RG functions owing to the UV finitcness of th e latter. 

The coordinate of the fixed point is determined by the condition that /?(<?*) = 0. From (3.1) and ( |3.2| ) we thus 
obtain: 



g* = (407r 2 e/3)(I + Ae) + 0(e 3 ), A = 2a 2 i/3<z 2 1 ~ -1.10. 



(3.3) 



The correction exponent lo is determined by the slope of the beta function at the fixed point, u> = (<?*)■ Thus from 
(3.2) and (3.3) we obtain the first and second terms of its e expansion: 



uj = 2e(l - Ae) + 0(e 3 ). 



(3.4) 



IV. TWO-LOOP CALCULATION OF THE KOLMOGOROV CONSTANT 

The Kolmogorov constant Ck can be defined as the (dimensionless) coefficient in the inertial-range expression 
S>x{f) = Cxitr) 2 ' 3 for the second-order structure function, predicted by the Kolmogorov-Obukhov theory and con- 
firmed by experiment . Here e is the mean energy dissipation rate and the n-th order (longitudinal, equal-time) 
structure function is defined as 

S n (r) = ([ip r (t,x + r)-ip r (t,x)] n ), ip r = ((pi-n)/r, r=\r\. (4.1) 

Using the exact relation S3 (r) = — 4er/5 that follows from the energy balance equation |l|,||, the constant Ck can be 
related to the inertial-range skewness factor: 

S = S a /S 3 2 /2 = -(4/5)C K 3/2 . 

Many studies have been devoted to the derivation of Ck within the framework of the RG approach; see Refs. 
[11-18]. In order to obtain Ck, it is necessary to augment t he s olution of the RG equation for S 2 by some formula 
that relates the amplitude Do in the random force correlator (2.4) to the physical parameter e. In particular, in jTl 12 



the first-order term of the e expansion for the pair correlator was combined with the so-called eddy damped quasi- 
Markovian approximation for the energy transfer funct ion , taken at e = 2. More elementary derivation, based on the 



exact relation between e and the function d,F{k) from (2^4) was given in |13J; see also [|3|]9[]. In spite of the reasonable 



agreement with the experiment, such derivations are not immaculate from the theoretical viewpoints. Their common 
flaw is that the relation between e and Dq is unambiguous only in the limit e — > 2: 

so that the coefficients of the corresponding e expansions appear in fact arbitrary; see the discussion in jl6| and 
Sec. 2.10 of ||. The ambiguity is a consequence of the fact that the notion of the Kolmogorov constant has no definite 
extension to the nonphysical range e < 2. 

The experience on the RG theory of critical behaviour shows that unambiguous e expansions can be written for 
universal quantities, such as critical exponents, normalized scaling functions and ratios of amplitudes in scaling laws 
The constant Ck extended to the range e < 2 as in |ll|-[l3[| involves a bare parameter, Dq, and hence is not 
universal. 

To circumvent this difficulty, we propose below an alternative derivation that relates Ck to an universal quantity 
and thus does not involve any relation between Dq and e, and calculate Ck to second order of the expansion in e 
(previous attempts have been confined with the first order). Consider the ratio 

Q(e) = rd r S 2 (r)/\S 3 (r)\ 2 / 3 = rd r S 2 (r)/(-S 3 (r)) 2 / 3 . (4.3) 

The operation rd r = rd/dr kills the constant contribution ((p 2 ) in S2 that diverges as A — » 00 for e < 3/2; in S3 such 
constant contributions are absent. 

Solving the RG equations for the quantities in Q(e) in the IR range (Ar 3> 1) for general < e < 2 gives 

S 3 (r) = D Q r- 3A -f 3 (e), rd r S 2 (r) = D^V 2 ^/ 3 f 2 (e), A v = 1 - 2e/3; (4.4) 
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see e.g. Thus the quantity Q{e) — f2/{—f?,) 2 ^ in ( [0| ) does not depend on D and can be calculated in the 

form of a regular e expansion. We calculated the scal ing functions /2,3 to the second order of the e expansion, which 



corresponds to the two-loop approximation in (p.3j), (3.4), and obtained 



Q(e) = (l/3)(20e) 1/3 [l + 0.525e + 0{e 2 )] . (4.5) 

It is worth noting that the s expansion for can be obtained not only from the direct perturbative calculation, but 
also from the exact expression 



that follows from the energy balance equation and in the limit e — > 2, along with the formula ( |4.2| ), reproduces the 
correct coefficient —4/5 (see above). 

The value of Q(e) at e = 2 determines the Kolmogorov constant and skewness factor through the exact relations 
C K = 6 • 1CT 2/3 Q(2), S = -[1.5 ■ Q(2)]- 3 / 2 , which follow from the definitions and the identity rd r r x = \r x for am 



A. Substituting the value of (4.5) we obtain Ck = 3.02, S = —0.15. If we retained only the first-order term in ( |4.5| ) 
we would have obtained Ck = 1-47 and S = —0.45. We also recall the experimental estimates recommended in [h 
C K w 1.9 and S w -0.28. 



V. CONCLUSION 



We have accomplished th e com plete two-loop calculation of the renormalization constant and RG functions for 
the stochastic problem (2.1)— (2.4) and derived the coordinate of the fixed point, the UV correction exponent u>, 
the Kolmogorov constant Ck and the inertial-range skewness factor S to the second order of the corresponding e 
expansions. The new point is not only the inclusion of the second-order correction, but also the derivation of Ck 
through an universal (in the sense of the theory of critical behaviour) quantity. 

Of course, one should have not expected that the second-order terms of the e expansions would be small in 
comparison to the first-order terms. The experience on the RG theory of critical behaviour shows that such corrections 
are not small for dynamical models (in contrast to static ones) and for amplitudes (in contrast to exponents); see ||. 
It is thus rather surprising that in our case the account of the two-loop contributions leads to reasonable changes in 
the results. 

Although the e 2 correction to uj in (3.4) is rather large, it does not change its sign and hence does not destroy the 
IR stability of the fixed point. 

The first-order approximation Ck — 1-47 underestimates, and the second-order approximation Ck = 3.02 overes- 
timates the conventional experimental value of the Kolmogorov constant Ck ~ 1.9 l]]- Thus the experimental value 
of Ck (and hence for S) lies in between of the two consecutive approximations. A similar situation is encountered 
for the well-known Heisenberg model B, where the analogue of the Kolmogorov constant is known exactly and lies 
between the first-order and second-order approximations given by the corresponding e expansion JTs| ] . If we assume, 
by the analogy with the Heisenberg model, that the (unknown) exact predictions for C'k and S lie between the first 
two approximations, we may conclude that our calculation has given a very good estimate for these quantities. 

The work was supported by the Nordic Grant for Network Cooperation with the Baltic Countries and Northwest 
Russia No. FIN-18/2001 and the GRACENAS Grant No. E00-3-24. 
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